Computational and Enzymatic Studies of Sartans in SARS-CoV-2 Spike RBD-ACE2 Binding: The Role of Tetrazole and Perspectives as Antihypertensive and COVID-19 Therapeutics

This study is an extension of current research into a novel class of synthetic antihypertensive drugs referred to as “bisartans”, which are bis-alkylated imidazole derivatives bearing two symmetric anionic biphenyltetrazoles. Research to date indicates that bisartans are superior to commercially available hypertension drugs, since the former undergo stronger docking to angiotensin-converting enzyme 2 (ACE2). ACE2 is the key receptor involved in SARS-CoV-2 entry, thus initiating COVID-19 infection and in regulating levels of vasoactive peptides such as angiotensin II and beneficial heptapeptides A(1-7) and Alamandine in the renin–angiotensin system (RAS). In previous studies using in vivo rabbit-iliac arterial models, we showed that Na+ or K+ salts of selected Bisartans initiate a potent dose–response inhibition of vasoconstriction. Furthermore, computational studies revealed that bisartans undergo stable binding to the vital interfacial region between ACE2 and the SARS-CoV-2 “receptor binding domain” (i.e., the viral RBD). Thus, bisartan homologs are expected to interfere with SARS-CoV-2 infection and/or suppress disease expression in humans. The primary goal of this study was to investigate the role of tetrazole in binding and the network of amino acids of SARS-CoV-2 Spike RBD-ACE2 complex involved in interactions with sartans. This study would, furthermore, allow the expansion of the synthetic space to create a diverse suite of new bisartans in conjunction with detailed computational and in vitro antiviral studies. A critical role for tetrazole was uncovered in this study, shedding light on the vital importance of this group in the binding of sartans and bisartans to the ACE2/Spike complex. The in silico data predicting an interaction of tetrazole-containing sartans with ACE2 were experimentally validated by the results of surface plasmon resonance (SPR) analyses performed with a recombinant human ACE2 protein.


Introduction
Since the advent of the Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2) pandemic in late 2019, intense research into the discovery and design of precision antiviral drugs was conducted [1][2][3][4][5][6][7][8][9]. The recent announcement of the Pfizer and Moderna antivirals, Paxlovid and Molnupiravir, respectively, are examples of ongoing research and development efforts in this area [10][11][12][13][14][15]. The drugs have different mechanisms of action, with Paxlovid inhibiting the main 3CL viral protease and Molnupiravir interfering with RNA-dependent RNA polymerase. Such antivirals are crucial armaments in fighting not only the current pandemic but also future outbreaks, since they are less expensive to manufacture and distribute than vaccine therapies, and they can be orally administered in a prophylactic regimen or post-infection to minimize severe symptoms and hospitalizations. Furthermore, the combination of drugs operating through different mechanisms may optimize the inhibitory effect through synergy [16]. Similar therapeutic approaches were applied to other viruses, such as the human immunodeficiency virus [17,18].
Our research team discovered and synthesized a new cohort of antihypertensive drugs that were predicted to stably dock at catalytic site and semi-stably at furin-cleavage site (residues 681-686) of the transmembrane angiotensin-converting enzyme 2 (ACE2), which is the primary receptor for attachment and entry of SARS-CoV-2 into susceptible host tissues [19][20][21][22][23][24]. We refer to this new class of drugs as "bisartans" since they share structural similarities to the "sartan" therapeutics (e.g., Losartan, Olmesartan, Telmisartan, Irbersartan, Valsartan, etc.) that currently dominate the global market for the treatment of hypertension and heart disease. Bisartans are bis-alkylated imidazole derivatives bearing dual symmetric anionic biphenyltetrazole moieties ( Figure 1) [22]. They exhibit a proclivity to interact with metallo-receptors/enzymes, particularly those harboring Zn 2+ cofactors. In general, the zinc-binding motif in metalloenzymes or metalloproteases is rich in histidines [25], which are excellent sources of pi-pi interactions with ligands containing phenyl groups such as sartans. Our computational modelling and in vivo experimental results (from rabbit-iliac arterial models) suggest that bisartans could be beneficial for the treatment of not only heart disease, diabetes, renal dysfunction, and related vascular illnesses, but also COVID-19. The in silico docking and molecular dynamics (MD) simulations revealed that several bisartan homologs (such as BisA, BisB, BisC, and BisD) [22] exhibited enhanced binding affinity for the ACE2/Spike protein complex (PDB ID: 6LZG) compared to other known (marketed) sartans, including the angiotensin receptor I blocker lisinopril (Figures 2 and 3) [22]. Bisartans undergo stable docking to the Zn 2+ -domain of the ACE2 catalytic site, as well as the critical interfacial region between ACE2 and the SARS-CoV-2 receptor-binding domain ( Figure 3). Additionally, bisartans exhibit semi-stable docking to the furin-cleavage site (residues 681-686) of the SARS-CoV-2 Spike protein required for viral entry into host cells. The docking studies suggest furin inhibition could occur by drug interaction with this arginine-rich site and, in particular, the P681R mutation which conveys enhanced infectivity to the Delta variant [19]. These findings strongly point to the potential of bisartans to function as a novel class of antiviral agents in addition to their antihypertensive action [21,22]. Int [22].   al. [22]. Docking  The major goal of our research was to further develop and experimentally evaluate in vitro and in vivo the efficacy of bisartans, not only as multifunctional SARS-CoV-2 antiviral agents but also as antihypertension and antibacterial drugs. In addition, we plan to design a range of bisartan structural analogs and/or pharmacophore modifications and computationally evaluate their behavior against additional human disease targets (receptors and enzymes), especially those containing metal cofactors. To achieve these objectives, we plan to focus on the power of the latest machine-learning chem/bioinformatics approaches to drug discovery and design, including the implementation of advanced algorithms offered by Molsoft, Schroedinger, and other commercial and open-source software developers. In addition, by the use of surface plasmon resonance technology, we aim to experimentally validate the binding capacity of sartans to the SARS-CoV-2 receptor ACE2.  [22]. Docking was performed using AutoDock VINA (YAMBER3 force field) with 100 runs for each ARB. The docking domain comprised a rectilinear cell with non-periodic (wall) boundaries 8 Å from any target atom. (A) Overview of 6LZG showing the major ACE2 channel where ARBs preferentially underwent docking. (B) The binding motif of bisartan BisA involved hydrophobic (green lines), pi-pi-resonance (red lines), and cation-pi (blue lines) interactions. Although the BisA conformational pose shown was proximate to the Zn 2+ cofactor and also exhibited the strongest binding (9.93 kcal/mol) compared to all other ARBs, it did not directly coordinate with the Zn 2+ cofactor. (C) The binding motif of BisB (8.25 kcal/mol) underwent coordination with Zn 2+ via pi-pi resonance (red lines) and cation-pi (blue lines) interactions. (D) VINA docking scores (bars) and dissociation constants (Kd; orange line) were calculated for the best poses of all 15 ARBs. With the exception of BisB, the bisartans (green-shaded bars) consistently exhibited stronger ACE2 binding compared to all other ARBs.

Results
The major goal of our research was to further develop and experimentally evaluate in vitro and in vivo the efficacy of bisartans, not only as multifunctional SARS-CoV-2 antiviral agents but also as antihypertension and antibacterial drugs. In addition, we plan to design a range of bisartan structural analogs and/or pharmacophore modifications and computationally evaluate their behavior against additional human disease targets (receptors and enzymes), especially those containing metal cofactors. To achieve these objectives, we plan to focus on the power of the latest machine-learning chem/bioinformatics approaches to drug discovery and design, including the implementation of advanced algorithms offered by Molsoft, Schroedinger, and other commercial and open-source software developers. In addition, by the use of surface plasmon resonance technology, we aim to experimentally validate the binding capacity of sartans to the SARS-CoV-2 receptor ACE2.

Computational Studies on the Interactions of Ligands with Spike Protein/ACE2
The effect of our generated ARBs and some well-known sartans on the SARS-CoV-2 spike protein/ACE2 region was investigated using docking and PLIP (PDB ID: 6LZG) [10.1016/j.cell.2020.03.045], and their binding affinity is shown at Table 1. The smaller molecules Co1 and Co2 failed to bind effectively and provided almost the same score, with the active site of the S-protein due to the limited number of interactions. This can be also excused due to the lack of aromatic rings that can provide the molecule to bind with the protein with π-stacking and π-cation interactions. Both molecules formed hydrogen bonds with 353A and 496B Gly, as well as hydrophobic bonds with 497B Phe and 505B Tyr.
In an effort to form interactions with aromatic rings, we decided to keep the trt groups in our molecules. Co3 and Co5 displayed a higher binding affinity of −9.4 and −9.71 kcal/mol, respectively. Compared to valsartan, olmesartan, eprosartan, and losartan (−8.53, −7.88, −6.98 and −7.74 kcal/mol), these molecules bound more effectively (Tables 1 and 2). Those molecules formed the biggest number of bonds compared to the other compounds due to their aromaticity. Although they did not form the same bonds with the S-protein, they formed similar ones. They both formed hydrophobic interactions with 505B Tyr, a π-stacking with 34A His, and a π-cation interaction with 403B Arg. Following the same pattern, we decided to generate Co7 and Co9 with the only difference from the previous group being the presence of a -CH 2 OH group, we thought that maybe more hydrogen bonds can be formed. However, that did not happen; on the contrary, the new compounds had a lower binding affinity (−6.8 and 6.81 kcal/mol) compared to Co3 and Co5. Again, all the compounds were aimed at the active site of S-protein; thus, all of them interacted with the same amino acids. Co7 and Co9 formed hydrogen bonds with 33A Asn and 393A Arg, hydrophobic interactions with 30A Asp, 33A Asn, 34A His, 37A Glu, 387A Ala, and 417B Lys. The interactions are shown in detail with their respective distances in Table 3.   Then, we removed the -trt group from the same molecules. As shown, by removing that group, the binding affinity of those molecules increased significantly in all the cases, with the most drastic increase in the case of Co9 and Co10, making the latter an excellent ligand (−11.4 kcal/mol). By removing the -trt from our previously best-performing molecule, Co5, the produced Co6 had an even lower binding energy (−12.1 kcal/mol). The resulting Co4 and Co8 also bound strongly to S-protein. By decreasing the size of the previous molecules, thus limiting their hydrophobic nature, we observed fewer hydrophobic bonding. However, it was also apparent that the number and quality of the forming hydrogen bonds increased. This was due to the fact that the -trt groups can limit the compounds' ability to be near several amino acids due to the steric effect, which, otherwise, could bind with them effectively. The interactions are shown in detail with their respective distances in Table 4. Table 4. The average values of binding affinity for the complexes of our generated compounds Co4, Co6, Co8, and Co10, (which do not contain the trt-group like their Co1-Co3, Co5, Co7, and Co9 counterparts) with COVID-19 S-Protein. The average values of binding affinity for the complexes of our generated compounds Co4, Co6, Co8, and Co10, (which do not contain the trt-group like their Co1-Co3, Co5, Co7, and Co9 counterparts) with COVID-19 S-Protein.

Co4
Co6 Co8 Co10 Hydrogen To elucidate the binding abilities of the new class of sartans to the ACE2/Spike protein complex, several binding experiments were performed for both the well-known sartans and the novel bisartan homologues. In the case of candesartan, eprosartan, irbesartan, and olmesartan ( Figure 4), the ligands bound to the catalytic center of ACE2 through strong pi-pi Interactions between the sartans' phenyl groups and histidines of the zinc-binding motif, including His 408 and 378. Additionally, salt bridges between carboxyl moiety and His 505 and Arg 273 of ACE2 stabilized the binding pose of eprosartan in the binding pocket. Tyr 515 of ACE2 played a critical role in the binding event of sartans due to its ability to interact via hydrogen bonds with nitrogen rings of ligands.

Enzymatic Studies Supporting the Interaction of Sartans with ACE2
To experimentally address the interaction of tetrazole-containing sartans with ACE2 protein, we performed SPR analysis using increasing amounts of the bisartan BV6 (containing two tetrazoles), the sartan losartan (containing one tetrazole and one hydroxyl) or its metabolite Exp3174 (losartan carboxylic acid containing one tetrazole and one carboxyl), and a recombinant human His-tagged ACE2 ( Figures 5-7). An efficient binding to ACE2 was shown for all sartans. The order of binding was BV6 > Losartan > Losartan carboxylic (Exp3174) and showed the superiority of tetrazolate compared to hydroxylate and carboxylate in binding to ACE2. For ACE2-Losartan complex formation, a kinetic constant was calculated with Biacore T200 Evaluation Software using the 1:1 binding model. The KD was 1.37

Enzymatic Studies Supporting the Interaction of Sartans with ACE2
To experimentally address the interaction of tetrazole-containing sartans with ACE2 protein, we performed SPR analysis using increasing amounts of the bisartan BV6 (containing two tetrazoles), the sartan losartan (containing one tetrazole and one hydroxyl) or its metabolite Exp3174 (losartan carboxylic acid containing one tetrazole and one carboxyl), and a recombinant human His-tagged ACE2 (Figures 5-7). An efficient binding to ACE2 was shown for all sartans. The order of binding was BV6 > Losartan > Losartan carboxylic (Exp3174) and showed the superiority of tetrazolate compared to hydroxylate and carboxylate in binding to ACE2. For ACE2-Losartan complex formation, a kinetic constant was calculated with Biacore T200 Evaluation Software using the 1:1 binding model. The KD was 1.37 × 10 −8 M (SE 1.5 × 10 −8 M).

Computational Approaches
Using state-of-the-art computational approaches, we identified a number of therapeutic targets for bisartans, namely three targets essential for viral infection and replication (i.e., ACE2, Furin, 3CLpro) [22]. These studies identified two additional therapeutic targets for bisartans, namely: (1) the zinc-dependent metallo-beta-lactamase from Bacillus cereus [20], and (2) Neprilysin, a zinc-dependent type II integral membrane peptidase belonging to the M13 family [20]. Beta-lactamases are responsible for imparting multiple antibiotic resistance to a range of human bacterial pathogens via hydrolysis of the beta-lactam ring of penicillin and its many derivative antibiotics, severely limiting therapeutic options. Neprilysin (NEP), on the other hand, is a regulatory enzyme that degrades beneficial vasoactive peptides, such as the atrial natriuretic peptide, the brain natriuretic peptide, bradykinin, adrenomedullin, and endothelin-15. Thus, the inhibition of NEP by bisartans is expected to lead to increased levels of these peptides. The combination of NEP suppression and angiotensin receptor inhibition (e.g., by bisartans) is known to be superior to the inhibition of either agent alone and leads to vasodilation and reduction in extracellular fluids via Na + ion excretion.

Hypertension and COVID-19 Mechanisms Are Similar
Cardiovascular disease is related to COVID-19 in terms of mechanisms that trigger infection [26][27][28][29]. The release of the cytokine storm in patients with severe pneumonia is related to the over-expression of toxic angiotensin II in the renin-angiotensin system (RAS) [30,31]. The role of RAS in autoimmune inflammation, in the pathogenesis of major human diseases, and as a therapeutic target, was extensively investigated [32][33][34][35]. The alteration from pro-inflammatory to anti-inflammatory secretion in autoimmunity was reported [36]. Clinical studies advocate in favor of the beneficial effects of ARBs in blocking infection [26,27,31]. Hypertension gained popularity among researchers due to its overrepresentation among COVID-19 patients [37]. Studies reported that hypertension is the most common co-morbidity observed in patients affected by SARS-CoV-2 [38][39][40]. The mechanisms that link pre-existing hypertension and COVID-19 are related to endothelial dysfunction and RAS imbalance. The conventional RAS (ACE/Ang II/AT1R) axis activation in parallel with nonconventional (ACE2/Ang 1-7/Mas) axis down-regulation was proposed to be the underlying factor leading to severe COVID-19 outcome in hypertension [41,42]. The imbalance that favors the pro-inflammatory state is proposed to be the center of COVID-19 pathophysiological mechanisms [43].

Anionic Groups of ANGII and ARBs Interact with Positive Sites of AT1R and ACE2
Recent crystal and in silico studies investigating the interactions between angiotensins and angiotensin receptor blockers (ARBs) with AT1R showed that sartans and bisartans which bear anionic tetrazolate and carboxylate groups as the pharmacophoric warheads bind to the 167Arg residue of AT1 receptor [28,29] and to the ACE2/RBD complex used by SARS-CoV-2 to infect cells [22]. In particular, olmesartan and carboxylic losartan (Exp3174) use the tetrazolate and carboxylate anionic groups to form the critical salt bridges with Arg 167 of AT1R and bisartans, bearing two tetrazolates, use both groups for the interaction with the receptor [19,21,22,44] (Figure 8). Aromatic interactions with Trp at position 84 strengthen the binding. As reported, the guanidino group of Ang II interacts with carboxyls of Asp363 and Asp282 of the AT1R and the anionic groups of the peptide are coordinated with Zn-bound ACE2. These studies showed the ability of bisartans to act as ARBs inhibiting vessel constriction in response to cumulative doses of ANGII and to reduce vasoconstriction due to the tight binding to Arg 167 of the AT1 receptor. The apparent targeting of multifunctional sites of the ACE2/RBD complex uncovered by in silico studies render sartans and bisartans bearing tetrazolate/carboxylate warheads promising drugs against COVID-19 [28,29,45,46].  [28,29]. Arg167 of the AT1 receptor interacts with tetrazolate and carbo ylate of olmesartan. (B) Interactions of Exp3174 (losartan carboxylic acid) and ACE2 in RBD spik ACE2 complex. Tetrazole of Exp3174 interacts with ACE2 Arg393.

Tetrazole of Sartans Interact with AT1R (Arg 167) and Spike 681-686 Arginines
Structure-activity, fluorescence, nuclear magnetic resonance, and computation studies led to a conformational model where the three aromatic amino acids (Tyr, H Phe) formed a tripartite ring cluster with the hydroxylate of tyrosine to trigger activi [47][48][49][50][51][52][53]. The model and the importance of the three aromatic amino acids for activity we supported by the design and synthesis of appropriate constrained cyclic analogues [5 using Barlos 2-chlorotrityl resin [55]. Alterations in the three critical aromatic residues angiotensin II deleted the agonist activity, confirming the importance of the aromatic re idues for the activity of the peptide. A charge relay system was suggested to operate f triggering the activity of AT1R analogous to the charge relay system of serine proteas [56]. Anionic interactions between ANGII and AT1R were in analogy to the interactio between ARBs and AT1R identified in the crystal structure of the ARB olmesartan wi  [28,29]. Arg167 of the AT1 receptor interacts with tetrazolate and carboxylate of olmesartan. (B) Interactions of Exp3174 (losartan carboxylic acid) and ACE2 in RBD spike-ACE2 complex. Tetrazole of Exp3174 interacts with ACE2 Arg393.

Tetrazole of Sartans Interact with AT1R (Arg 167) and Spike 681-686 Arginines
Structure-activity, fluorescence, nuclear magnetic resonance, and computational studies led to a conformational model where the three aromatic amino acids (Tyr, His, Phe) formed a tripartite ring cluster with the hydroxylate of tyrosine to trigger activity [47][48][49][50][51][52][53]. The model and the importance of the three aromatic amino acids for activity were supported by the design and synthesis of appropriate constrained cyclic analogues [54] using Barlos 2-chlorotrityl resin [55]. Alterations in the three critical aromatic residues of angiotensin II deleted the agonist activity, confirming the importance of the aromatic residues for the activity of the peptide. A charge relay system was suggested to operate for triggering the activity of AT1R analogous to the charge relay system of serine proteases [56]. Anionic interactions between ANGII and AT1R were in analogy to the interactions between ARBs and AT1R identified in the crystal structure of the ARB olmesartan with AT1R [28,29]. The crystal structure shed light on the interactions and the critical role of the AT1R Arg167 bound to the two anionic pharmacophoric tetrazole and carboxylate groups of olmesartan. In silico studies showed that the ionic groups of sartans and bisartans, in particular tetrazole, bound to the positive Arg167 guanidino group of AT1R and the arginines of the ACE2/RBD complex [22]. These studies, furthermore, showed interactions of sartan tetrazole with arginines in the rich Arg cleavage site 681-686, thus preventing cleavage and infectivity [19][20][21][22]. Figure 9 shows the structures of designed and synthesized imidazole-based ligands, protected by trityl group and unprotected, for studying their interaction with the RBD/ACE2 complex. Computational studies showed bisartan bisA (Co4) to be best docked with the RBD/ACE2 complex [22]. The free tetrazole of the investigated bisartans (Co4, Co6, Co8, Co10) bound stronger compared to trityl protected tetrazoles (Co3, Co5, Co7, Co9).

Imidazole and Benzimidazole Based Sartans
There are two subclasses of sartans, imidazole and benzimidazole-based sartans. Examples of imidazole scaffold sartans are losartan, Exp3174 (carboxylic losartan), olmesartan, eprosartan, and all those sartans built on the imidazole ring. Examples of benzimidazole scaffold sartans are candesartan, telmisartan, Azilsartan, and those sartans built on benzimidazole moiety which fuses histidine and phenyl rings. Telmisartan was found to protect hypertensive patients infected by SARS-CoV-2 [57,58] and in general, the clinical trials advocate in favor of the beneficial effects of ARBs in blocking infection [26,27,31,59]. These studies showed that morbidity and mortality rate was lower in hypertensive patients infected by SARS-CoV-2 who were prescribed RAS and, in particular, ARBs inhibitors when compared to patients not taking these drugs [57,[60][61][62]. Hypertensive patients not taking RAS antihypertensives are more vulnerable to developing serious complica-Scheme 1. The scheme depicts the synthesis of imidazole-based bisartans bearing two symmetric biphenyl tetrazoles and butyl/CH 2 OH groups at positions 2 and 4, respectively. Figure 10 shows the crystallographic grids of the generated ligands (Co1-Co10) as well as the tested sartans valsartan, olmesartan, eprosartan, and losartan. The latter are depicted at Figure 11. Figure 9 shows the structures of designed and synthesized imidazole-based liga protected by trityl group and unprotected, for studying their interaction with RBD/ACE2 complex. Computational studies showed bisartan bisA (Co4) to be docked with the RBD/ACE2 complex [22]. The free tetrazole of the investigated bisar (Co4, Co6, Co8, Co10) bound stronger compared to trityl protected tetrazoles (Co3, Co7, Co9). Figure 9. Structures of the used ligands. Bisartans are shown with protection of tetrazole with group (Co3, Co5, Co7, Co9) and deprotected (Co4, Co6, Co8, Co10) (Scheme 1 shows the synt steps). Figure 10 shows the crystallographic grids of the generated ligands (Co1-Co1 well as the tested sartans valsartan, olmesartan, eprosartan, and losartan. The latte depicted at Figure 11.  (Co3, Co5, Co7, Co9) and deprotected (Co4, Co6, Co8, Co10) (Scheme 1 shows the synthesis steps).

Imidazole and Benzimidazole Based Sartans
There are two subclasses of sartans, imidazole and benzimidazole-based sartans. Examples of imidazole scaffold sartans are losartan, Exp3174 (carboxylic losartan), olmesartan, eprosartan, and all those sartans built on the imidazole ring. Examples of benzimidazole scaffold sartans are candesartan, telmisartan, Azilsartan, and those sartans built on benzimidazole moiety which fuses histidine and phenyl rings. Telmisartan was found to protect hypertensive patients infected by SARS-CoV-2 [57,58] and in general, the clinical trials advocate in favor of the beneficial effects of ARBs in blocking infection [26,27,31,59]. These studies showed that morbidity and mortality rate was lower in hypertensive patients infected by SARS-CoV-2 who were prescribed RAS and, in particular, ARBs inhibitors when compared to patients not taking these drugs [57,[60][61][62]. Hypertensive patients not taking RAS antihypertensives are more vulnerable to developing serious complications of COVID-19. The imbalance in the renin-angiotensin system (excess of toxic angiotensin II against beneficial heptapeptides A(1-7), Alamantine) is responsible for hypertension and COVID-19. ARBs upregulate ACE2 which upgrades angiotensin II are proven to be beneficial against SARS-CoV-2 [63,64]. Scheme 1 shows the synthetic scheme for bisartan imidazole-based ligands Co4, Co6, Co8, and Co10. Int

Imidazole and Benzimidazole Based Sartans
There are two subclasses of sartans, imidazole and benzimidazole-based sartans. Examples of imidazole scaffold sartans are losartan, Exp3174 (carboxylic losartan), olmesartan, eprosartan, and all those sartans built on the imidazole ring. Examples of benzimidazole scaffold sartans are candesartan, telmisartan, Azilsartan, and those sartans built on benzimidazole moiety which fuses histidine and phenyl rings. Telmisartan was found to protect hypertensive patients infected by SARS-CoV-2 [57,58] and in general, the clinical trials advocate in favor of the beneficial effects of ARBs in blocking infection [26,27,31,59]. These studies showed that morbidity and mortality rate was lower in hypertensive patients infected by SARS-CoV-2 who were prescribed RAS and, in particular, ARBs inhibitors when compared to patients not taking these drugs [57,[60][61][62]. Hypertensive patients not taking RAS antihypertensives are more vulnerable to developing serious complications of COVID-19. The imbalance in the renin-angiotensin system (excess of toxic angiotensin II against beneficial heptapeptides A (1-7), Alamantine) is responsible for hypertension and COVID-19. ARBs upregulate ACE2 which upgrades angiotensin II are proven to

Possible Target of Sartans as Revealed by Enzymatic Assays
In the current study, we selected to study the SPR binding to the human ACE2 protein of free sartans, not trityl protected, which were predicted to bind better to the ACE2-Spike RBD complex in computational studies, as depicted in Table 1 (best scoring pose bound to ACE2-Spike RBD complex). In particular, we used the bisartan BV6 (containing two tetrazoles), the sartan losartan (containing one tetrazole and one hydroxyl), and losartan carboxylic acid, Exp3174 (containing one tetrazole and one carboxyl), for the ACE2 study. The results presented showed an efficient binding of the three sartans. The order of binding BV6 > Losartan > Losartan carboxylic (Exp3174) suggests the superiority of tetrazolate compared to hydroxylate and carboxylate in binding to ACE2 in line with our previous studies, where BV6 was the best binder in computational studies (Ridgway et al. [22]). On the other hand, both sartans were found to inhibit the enzymatic activity of a recombinant soluble ACE2 only slightly (less than 10%) at a concentration of 100 micromolar, which suggests that sartans bind to amino acids at an area remote from the catalytic ACE2 Zn 2+ center, not directly to the active site. Taken together, the binding and enzymatic studies suggest that the specific sartans do not selectively target the active site of the ACE2 enzyme but rather its interaction domain with the RBD of the viral S1 protein.

Clinical Perspectives of Bisartans as COVID-19 Antivirals
The recent announcement of the Pfizer and Moderna antivirals, Paxlovid and Molnupiravir, triggered ongoing research for new antivirals against COVID-19 based on the action of their mechanism [11][12][13][14][15][16]. These drugs use different mechanisms of action to exert their antiviral activities against SARS-CoV-2 spike protein, first by blocking 3CLpro and second by interfering with RNA-dependent RNA polymerase, resulting in transcription error accumulation. Paxlovid's active ingredient nirmatrelvir (PF-07321332) through the warhead nitrile group forming a covalent bond with thiol of the catalytic dyad Cys145-His41 of 3CLpro inactivates the catalytic site and blocks cleavage of the spike protein, thus preventing infection [10,14,65]. Our in silico studies showed that bisartans block the spike cleavage as arginine blockers through the three cell entries (ACE2, furin, 3CLpro) rendering them promising drugs for treating COVID-19 [21,22].

In Silico Methodology and Ligand Preparation
All compounds were sketched in the Chem3D 15.0 module of ChemOffice 15.0 and converted into SMILES [22]. This was accomplished using Maestro MacroModel 10.8 (Schroedinger, LLC, New York, NY, USA) where all hydrogens were added and molecules were subjected to complete structure minimization. For the minimization, water was chosen as a solvent, OPLS3 [66] force field, and the algorithm Polak-Ribiere (PRCG, convergence value 0.01 kcal mol −1 Å −1 ) were used. The pH of the in silico studies was set to 7.4 based on the pH of human blood, which made our environment slightly basic. As a result, the protons of the tetrazole rings as well as the protons of the carboxylic acids are not present. The resulting files were saved in pdb format [67]. In AutoDock 4.0, flexible torsions were assigned allowing 100 conformations [68,69] and the acyclic dihedral angles were allowed to rotate freely. The files were then saved in the pdbqt file format for further analysis [67].

Molecular Docking
Molecular simulation studies of our compounds into the protein targets were carried out using the open-source program Autodock 4.0 included in Auto-Dock Tools 1.5.6. For this study, the crystal structures of the proteins were extracted by the Protein Data Bank. These results were compared to candesartan, a known ARB, from our previous study, as well as valsartan, olmesartan, and eprosartan [63]. The X-ray crystal structures of the spike receptor domain complexed with ACE2 (PDB ID: 6LZG) were downloaded from the RCSB PDB (Protein Data Bank) database [70,71] [72][73][74].
The best-docked poses, with both lower binding energies and stronger interaction patterns, were derived from the docking results and were visualized with PyMOL and the protein-ligand interaction profiler (PLIP; retrieved from https://pymol.sourceforge. net/overview/index.htm (accessed on 30 October 2022), PyMOL version 2.0, Schrödinger, Inc., New York, NY, USA) [75]. The PLIP was also used to determine the interactions between the ligands and the proteins [76]. The X-ray crystal structures of the selected SARS-CoV-2 proteins were downloaded from Protein Data Bank as pdb files, and Open Babel software was used to convert the protein files into pdbqt types. The co-crystallized ligands were removed; thus, 2-acetamido-2-deoxy-beta-D-glucopyranose of spike protein (PDB ID: 6LZG) [77]. The protein structures were refined for heteroatoms and water molecules to demarcate the active sites of the proteins. The hydrogen atoms and the nonpolar hydrogens were merged, and Gastgeiger and Kollman charges were added. In each receptor, the grid was set around its active site for site-specific docking to be performed [78]. The Lamarckian genetic algorithm (GA), which uses the AMBER force field to run the docking between the receptor and the ligand, was used for docking in combination with the grid-based energy evaluation method with the default parameters (GA: 5,000,000-energy evaluations and 175-population size). The program was run for a total number of 100 genetic algorithm runs. When using Autodock 4.0, the stability of the ligand/protein complex, which shows how efficiently can the ligand binds to the protein, is depicted with low energy. The lower the energy, the more efficient the ligand binds, which is depicted by the increased number of hydrogen bonds, non-covalent Van Der Waals, and hydrophobic interactions.

Docking Parameters
For in silico SARS-CoV-2 protein experiments with the selected ligands, the active site of each receptor was targeted. For the complex spike protein, with ACE2 as its predominant receptor (PDB ID: 6LZG) [77], the grid box was set with a spacing of 0.420 Å and dimensions of 77 Å × 80 Å × 100 Å, centering around residues mentioned above. In silico molecular docking analysis was performed with Autodock 4.0 [67]. The ligands produced pdbqt files, which were converted to pdb, using Open Babel to be compatible with the visualization platforms. Then, the pdb files were visualized and analyzed by PyMOL and PLIP [75,79].

Organic Synthesis of Bisartans
Starting materials were purchased by Aldrich (Patras, Greece) and were used as received. The 1H-NMR and 13C-NMR spectra were recorded on a Bruker Avance DPX spectrometer at 400.13 MHz and 161.76 MHz, respectively. Chemical shifts are given in δ values (ppm) using tetramethylsilane as the internal standard and coupling constants (J) are given in Hertz (Hz). HPLC analysis was performed on an Alliance Waters 2695 equipped with a Waters 2996 Photodiode Array Detector UV-Vis, using the XBridge Waters C18 column (4.6 × 150 mm, 3.5 µm) as stationary phase and a gradient of H 2 O/CH 3 CN, both containing 0.08% TFA as mobile phase. Electrospray-ionization mass spectra (ESI-MS) were obtained on a UPLC (ultra-performance liquid chromatography) equipped with SQ detector AcquityTM by Waters. Analytical TLC was performed on silica gel 60 F254 plates (Merck, Darmstadt, Germany) and visualized by UV irradiation. Silica gel 60N (particle size 0.04-0.063 mm) was used for flash column chromatography.
General procedure for the bis-alkylations: A solution of 4(5)-butylimidazole (1) or 2-butylimidazole (2) (0.8 mmol, 0.1 g) and 5-(4 -(bromomethyl)biphenyl-2-yl)-2-trityl-2Htetrazole (1.69 mmol, 0,94 g) in 50 mL of dichloromethane (DCM) was refluxed for 24 h and the reaction was monitored by RP-HPLC (70% CH 3 CN in H 2 O to 100% CH 3 CN, in 30 min). The reaction mixture was diluted with 100 mL of dichloromethane and the organic solution was washed with 50 mL aq. KOH 1N, water (2 × 50 mL), and brine (50 mL). The organic layer was dried with over anhydrous Na 2 SO 4 , filtered and dichloromethane was evaporated in vacuo. In the oily residue, ethyl acetate (20 mL) and then, diethylether (100 mL) were added and a solid was precipitated. After filtration, the desired product (3) or (4) was afforded, in high yield and enough purity, in order to be used as it was for the next steps. Data for compound (3) 3 CN, in 30 min) and the mixture was diluted with chloroform (200 mL). The organic layer was separated and washed with 10% w/v aq. citric acid (50 mL), water (50 mL), and brine (50 mL), dried with over anhydrous Na 2 SO 4 , filtered and chloroform was removed by vacuum evaporation. Diethylether (100 mL) was added to the oily residue and a white solid was precipitated, and, after filtration, the desired product (5) or (6) was afforded, in high yield and purity, and was used as it is to the deprotection reaction. Data for compound (5) General procedure for the triphenylmethyl group removal: A quantity of a tetrazole Trt-protected bis-alkylated solid compound was dissolved in a 50% v/v trifluoroacetic acid (TFA) in dichloromethane solution and a few drops of triethylsilane (TES) were added, until decoloration of the mixture was observed. The resulting solution was stirred for 1 h at ambient temperature and then, was evaporated in vacuo. After diethylether addition, a solid was precipitated, which was obtained by filtration as an amorphous solid, in the form of the TFA salt of the respective protected compound. Data for compound (7): By deprotection of 0.2 g (0.17 mmol) of compound (3): White amorphous solid, 0.14 g (90% yield). RP-HPLC (20% CH 3

ACE2 Protein Purification
The human ACE2 protein, tagged with His-tag and Strep-tag, was overproduced in Sf 9 insect cells. For purification, 107 cells were used. The ACE2 tagged protein was purified with the usage of Strep-Tactin Superflow Plus resin (Qiagen, Hilden, Germany), according to the manufacturer's protocol. Cells were lysed in the NP buffer (50 mM NaH2PO4 pH 8, 300 mM NaCl, 10% (v/v) glycerol, 1% (v/v) Triton X-100) containing proteases inhibitors (Pierce™ Protease Inhibitor Mini Tablets). The clear cellular lysate was incubated with the resin for 1 h with rotation. Next, the resin was washed with NP buffer and ACE2 protein was eluted with NP buffer containing 10 mM biotin hydrazide.

Surface Plasmon Resonance Analysis
SPR experiments were performed with Biacore T200, (Cytiva, Malborough, MA, USA) equipment [80]. The purified ACE2-tagged protein was covalently immobilized on a CM5 Sensor Chip (Cytiva, Malborough, MA, USA) via the primary amine group, in a 10 mM sodium acetate buffer, pH 4.5, as described by Zhu and co-workers. The protein was immobilized to a response level of 7700 RU. Experiments were run at 25 • C, and the running buffer was PBS-P (Cytiva, Malborough, MA, USA) with 2% DMSO. The sartans were prepared running buffer in a series of concentrations from 0.6 to 180 µM and flowed over the immobilized ACE2 protein at a flow rate of 30 µL/min, contact time of 120 s, and dissociation time of 300 s. The sensor chip surface was regenerated using 1.5 M NaCl and 10 mM glycine, pH 1.5. After each analysis, an additional wash with 50% DMSO solution was performed. The results are presented as sensorgrams obtained after subtracting the background response signal from a reference flow cell and a control experiment with buffer injection. For ACE2-Losartan interaction, a kinetic constant was calculated with Biacore T200 Evaluation Software using the 1:1 binding model.

Conclusions
The present in silico study revealed the significant role of the tetrazole of sartans and bisartans for binding to the ACE2/Spike complex. The protection of warhead tetrazole with trityl group reduced the affinity of sartans to the complex, confirming the importance of negative tetrazolate for binding. SPR binding analysis confirmed the interaction of tetrazole-containing sartans with ACE2. Imidazole and benzimidazole scaffold-based sartans with tetrazole and carboxylate, as warheads, are promising compounds for treating COVID-19 and are worthy of further investigation.